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TECHNICAL  NOTE  288? 

ON  THE  STABILITr  OF  THE  LAMINAR  MIXING  REGION 
BETWEEN  TWO  PARALLEL  STREAMS  IN  A GAS 
By  C.  C.  Lin 

SUMARY 


The  stability  of  the  m-tyt-ng  of  tvro  streams  vas  studied  hoth  for 
the  Interest  In  the  problem  and  for  clarifying  certain  points  In  the 
basic  stability  theory.  It  is  shown  that,  when  the  relative  speed  of 
the  two  parallel  streams  exceeds  the  sum  of  their  velocities  of  sound, 
BTibsonic  oscillations  cannot  occnir  and  the  mixing  region  may  be  expected 
to  be  stable  with  respect  to  small  disturbances.  It  is  further  shown 
that,  when  viscosity  and  heat  conductivity  are  neglected,  if  the  flow 
nan  execute  a small  neutral  stibsonic  oscillation  of  finite,  wave  length, 
it  can  also  execute  self -excited  oscillations  of  longer  wave  lengths 
fl.rif=i  damped  osciHations  of  shorter  wave  lengths. 

Rigorous  developments  of  the  mathematical  theory  of  asymptotic 
solutions  confirm  previoiis  methods  of  solution  of  the  stability  equa- 
tions in  a conpressible  fliold.  This  theory  also  shows  that,  at  hl^ 
Reynolds  numbers,  the  damped  oscillations  in  a strictly  parallel  main 
flow  have  a structure  similar  to  that  of  the  vorticity  field  in  f\iUy 
developed  turbulent  flow. 

Sample  calculations  are  also  included  exhibiting  various  quantita- 
tive properties  of  these  small  oscillations. 


INTRODUCTION 


The  miYing  of  two  parallel  streams  of  gas  occurs  in  a number  of 
cases.  An  interesting  example  is  furnished  by  the  slip  stream  in  a 
three-shock  configuration.  It  has  long  been  suggested  that  such  laminar 
TTiiviTig  zones  could,  at  sufficiently  high  speeds,  be  stable  with  respect 
to  small  disturbances  although  they  are  known  to  be  very  unstable  at  low 
speeds.  The  purpose  of  the  present  investigation  is  to  find  out  some 
of  these  stability  characteristics.  Apart  from  the  development  of  the 
general  theory,  there  are  included  the  calculations  of  the  neutral  and 
unstable  oscillations,  the  extent  of  the  amplification,  and  other  related 
properties . 
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The  "basic  equations  for  the  stvuiy  of  small  disturhances  in  the 

houndary  layer  of  gases  have  heen  given  in  reference  1.  However, 
the  development  of  the  theory  there  has  specific  reference  to  the  case 
of  a layer  near  a solid  houndary.  As  pointed  out  in  references  2 and  3^ 
the  stahillty  theory  for  a mlylng  zone  in  an  incon5>ressihle  fluid  dif- 
fers from  that  for  a layer  near  a solid  houndary  in  that  solutions  of 
the  exponentlEil  type  Eire  uniraportant  - This  leads  to  the  concltision  that 
the  effects  of  viscosity  and  heat  conduction  are  negligible,  except  at 
very  low  Reynolds  nvmibers,  in  determining  the  characteristics  of  the 
oscillations.  To’  confirm  this  point,  a rigorous  mathematical  theory  of 
asyn5)totlc  solutions  was  developed  for  the  cong)ressihle  case  similar  to 
that  indicated  in  reference  4 for  the  incongpressihle  case.  In  view  of 
the  mathematical  interest  involved,  it  was  decided  that  this  basic  part 
of  the  present  investigation  would  he  published  separately  in  mathematiCEil 
Journals  (see  references  5 and  6),  and  only  the  main  results  and  their 
physlcEil  significance  will  he  presented  here. 

As  in  the  case  of  an  inconpressihle  fluid,  the  "invlscid"  case  is 
eiqjected  to  he  chEiracterlstic  of  the  behavior  of  the  disturhances  at 
moderately  large  Reynolds  numbers.  Ifost  of  the  studies  are,  therefore, 
made  in  the  invlscid  case.  However,  the  interpretation  of  the  Inviscld 
case  must  he  subjected  to  the  same  care  as  in  the  incompressible  case} 
that  is,  in  the  case  of  damped  disturhances,  the  differential  equation 
of  the  invlscid  flow  may  not  he  regsirded  as  valid  throughout  the  real 
eixls.  There  is  a finite  vLscoris  region  even  in  the  limit  of  vanishing 
viscosity.  The  complex  conjugate  of  the  amplified  disturbance  is  cer- 
tainly a solution  of  the  invlscid  equation,  hut  it  is  not  a limiting 
solution  of  the  complete  viscous  equation.  This  behavior  of  the  invlscid 
solution  reminds  one  of  the  vorticity  structure  of  fully  developed  tur- 
bulent flow  as  found  by  Batchelor  Eind  Townsend  (reference  T)  • 

The  nonexistence  of  subsonic-  disturhances  is  usually  Eissoclated 
with  the  stability  of  the  parallel  flow.  There  seems  to  he  some  basis 
for  doing  this,  eilthough  the  role  of  supersonic  disturhances  has  never 
heen  fully  clarified.  It  is  easy  to  see  that,  for  certain  combinations 
of  the  properties  of  the  two  streams,  it  is  impossible  to  have  a sub- 
sonic disturbance  relative  to  both.  Under  such  conditions,  one  may 
expect  stability.  These  conditions  for  stability  are  developed  herein 
and  are  expected  to  hold.  Irrespective  of  the  -viscous  effects. 

Applying  the  theory  of  stability  in  the  invlscid  case,  one  can 
further  nsirrow  down  the  possible  range  of  instability.  This  will  depend 
upon  the  velocl-by  and  temperature  distributions  in  the  sheer  zone.  In 
the  present  work,  calculations  are  made  for  gases  with  Preindtl  number 
equal  to  unity.  Although  the  condition  of  equal  total  enthalpy  in  the 
two  streams  is  also  tised,  it  is  shown  that  this  restriction  cem,  be 
immediately  removed  by  considering  a moving  observer.  It  is  found 
that  the  condition  of  stability  thus  found  does  not  differ  very  much 
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from  that  found  above  from  general  considerations.  Thus,  it  may  he 
surmised  that  the  exact  distributions  of  temperature  and  velocity  have 
only  a secondary  influence  on  the  stability  characteristics  in  the 
mixing  zone.  Thus,  the  approximations  used  in  the  present  calculations 
of  the  baaic  velocity  and  tenrperature  distributions  CEinnot  influence 
the  final  results  to  any  appreciable  extent. 

The  neutral  disturbances  are  of  two  kinds:  (l)  A steady  deviation 

and  (2)  an  oscillation  of  finite  wave  length.  Thus,  there  are  two 
branches  of  the  neutral  curve  at  infinite  Reynolds  mmibers.  They  may 
be  expected  to  join  together  at  low  Reynolds  numbers  enclosing  a region 
of  instability.  Calculations  of  neutral  and  amplified  oscillations  are 
carried  out  in  a nimiber  of  cases  with  one  stream  at  rest. 

This  investigation,  carried  out  at  the  ifessachusetts  Institute  of 
Technology,  was  sponsored  by  and  conducted  with  the  financial  assistance 
of  the  Rational  Advisory  Committee  for  Aeronautics.  The  author  is 
Indebted  to  Mr.  D.  W.  Dunn  for  his  valuable  suggestions  and  help  in  the 
preparation  of  the  final  version  of  the  report  and  to  Miss  Diana  Mason 
and  Mr.  W.  V.  Caldwell  for  their  help  in  making  the  numerical  calculations. 


STEADY  FLOW  IN  THE  LAMINAR  LAYER  BETWEEN 
TWO  PARALLEL  STREAJB 


The  basic  steady  flow  mder  discussion  is  a boundary- layer  flow  with 
no  body  forces  and  no  pressure  gradient.  The  basic  equations  are  (see 
list  of  symbols  in  the  appendix) : 
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the  pressiire  helng  a constant  throughout  the  field.  In  the  case  of  a 
homogeneoTis  Incompresslhle  fluid,  with  the  streams  at  the  same  tempera- 
ture, the  ten5)erature  may  he  tahen  as  constant  throughout  the  field, 
and  the  first  two  equations  can  he  solved  for  the  velocity  distribu- 
tions u*(x*,y*)  and  v*(x*,y*) . In  the  case  of  a conipressihle  gets, 
the  integration  has  to  he  carried  out  for  individual  cases.  However, 
if  the  Pimidtl  number  Cpp*/k*  is  equal  to  unity,  it  is  known  that 

there  is  a quadratic  relation  of  the  type 


CpT*  + i u*^  = A + Bu* 


ik) 


between  the  temperature  and  the  velocity,  and  one  is  again  essentially 
dealing  with  two  distributions  u*(x*,y*)  and  v*(x*,y*).  Indeed,  the 
constants  A and  B are  given  in  terms  of  the  conditions  in  the, 
parallel  streams  as  follows; 


B ' 1(^2  + + <=p(T2  - Ti)/(02  - %) 

A = - I + <=p(TiD2  - ’^2“l)/(U2  " %) 


(5) 


If  the  total  enthalpy  in  the  two  streams  1 and  2 is  the  same,  that  is, 
if 


CpTi  + I = CpTg  + i V/ 


(6) 


then  B = 0,  and  the  total  enthalpy  is  constant  throughout  the  whole  field; 

CpT^  + i u*^  = C 


(7) 
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OOils  is  a particularly  simple  case,  to  whlcti,  however,  all  other  cases 
can  he  easily  reduced.  This  is  done  by  rewriting  equation  (!<-)  into 
the  form 


CpT*  + i(u*  - B)2  = A + b2/2 
2 


(8) 


One  need  only  consider  an  observer  moving  with  the  speed  B and  con- 
sider the  relative  velocity  u*  - B.  0?hus,  the  solution  for  general 
specified  values  of  U2_,T-|_,  and  Ug,T2  can  be  derived  fix)m  the  iso- 

energetic  solution  with  boundary  conditions  U2_  - B,Tjl  and  Ug  - B,Tg 
by  simply  adding  the  constant  B to  the  u*-component  of  the  velocity.^ 

In  the  following  disciissions,  isoenergetic  basic  solutions,  will  be 
referred  to  oftenj  however,  it  should  be  kept  in  mind  that  by  the  con- 
sideration of  a moving  observer  the  general  case  may  be  obtained.’  This 
transformation  is  not  restricted  to  the  steady  flow  but  applies  to  the 
consideration  of  the  dlstinr'bances  as  well.  0?hus,  if  aU  the  cases  of 
constant  total  enthalpy  are  calculated,  all  the  other  cases  are  also 
known. 

So  far,  the  viscosity  coefficient  may  depend  on  the  absolute  tem- 
perature in  any  manner.  If  there  is  direct  proportionality  of  these 
q.uantitles,  the  solution  in  the  conrpressible  case  can  be  expressed  in 
terms  of  that  in  the  incompressible  case.  OJhese  relations  are  well- 
kno-v/n  and,  in  the  following  discussion,  only  the  results  relevant  to 
this  case  will  be  given. 


Incompressible  Case 
For  the  incompressible  case: 


^ ^Vi(UiX*(Tif'  - f) 


(9) 


lit  is  to  be  noticed  that  in  certain  cases  with  T2/Tj_  < 1 the 
boundary  value  Ug  - B for  the  corresponding  isoenergetic  problem  may 
be  negative  while  Up  - B is  still  positive.  Thus  the  isoenergetic 
problem  may  not  be  physically  significant.  However,  for  the  purposes 
of  theoretical  analysis,  this  point  is  not  important. 
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where 


and"  f(T])  satisfies  the  differential  eqmtion 


ff " + 2f " ’ = 0 


(11) 


with  the  conditions 


f ’(ti)  — > 1 as  Tj  — > » 


f'(Tl)  >5^  as  T] 


(12) 


A third  condition  is  arbitrary  up  to  a translation  along  the  T)-axls . 

For  a typical  scale,  the  momentxmi-houndary- layer  thickness  9*  may 
he  introduced,  which  is  defined  hy 


/OD 

- u*^(u*  - Ug)  dy* 


(13) 


Then 


0 = 


f’(Tl) 


U2 


•U-, 


dll 


(llv) 


for  an  inconpressihle  fluid.  It  will  he  seen  later  that  the  same  for- 
mulas apply  for  a corapressihle  fluid  in  the  isoenergetic  case. 

Calculations  made  hy  Gortler  (reference  8)  for  the  tirrhulent  ml ng 
region  can  he  easily  adopted  for  the  pvupose  at  hand.  There  is  only  a 
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sllglit  difference  in  the  method  of  representation.  To  convert  his 
function  F'(|)  to  the  notation  of  this  report,  the  following  rela- 
tions shovild  he  used 


f«(n)  = 

1 + X 

Ti  = 2\Jl  + X I 

where 


Ui  +U2 

The  converted  resiilts  are  given  in  table  I and  figure  1. 


(15a) 


(15b) 


Compressible  Case 

For  the  congoressible  ca^e  with  constant  total  enthalpy,  p*  a T*, 

where 


u*  = U3^f'(^) 


(16) 


(17) 


(18) 


where 


T]  = 


8 


MCA  IN  2887 


The  momentum  thickness  6*  can  he  defined  hy  the  same  equation  (equa- 
tion (13))  and  can  he  shown  to  have  the  same  value  , (eqmtion  in 

dimensionless  form.  Velocity  and  temperature  dlstrihutions  in  the  case 
U2  = 0 have  heen  calculated  for  several  Mach  nuaibersj  these  are  tabu- 
lated in  table  H and  plotted  in  figure  2. 

A GENERAL  RESTRICTION  FOR  EXISTENCE  OF 
NEUTRAL  SUBSONIC  DISTURBANCES 


Before  going  into  the  general  theory  of  stability,  a preliminary 
discussion  of  the  stability  of  the  mixing  zone  will  be  given.  It  will 
be  found  below  that  the  results  thus  obtained  give  quite  an  adequate 
description  of  the  general  dependence  of  the  stability  of  the  mixing 
zone  on  the  Ifech  numbers  of  the  streams.  It  will  be  shown  that,  if 
the  average  Mach  number  It  of  the  relative  motion  (defined  by  equa- 
tions (25)  and  (26))  exceeds  the  value  2,  subsonic  disturbances  in  the 
sense  of  reference  1 cannot  exist,  and  the  mixing  zone  may  be  expected 
to  be  stable. 

Consider  two  parallel  streams  at  speeds  and  Ug,  tempera- 

tiores  T-]_  and  T2,  and  Mach  numbers  and  M2-  For  definiteness, 

tsike  U2_  - U2  > 0.  Consider  an  observer  moving  with  the  speed  U2. 
Then  the  streams  appear  to  have  speeds  Uq  - U2  and  0,  while  their 
teraperat-ures  are  obviously  not  changed.  If 


c'  = c - U2 


(19) 


denotes  the  speed  (relative  to  the  moving  observer)  of  a wavy  motion 
propagating  in  the  direction  of  the  stream  (c'  positive  or  negative), 
then  the  conditions  for  subsonic  distvirbances  are 


l(Ul  - Uj)  - c'l  < 


(20) 


jc  I < ag  (21) 

where  a^^  »tiH  are  the  velocities  of  sound  in  the  streams.  It  is 

clear  that  subsonic  disturbances  can  always  exist  if  - U2  is  less 
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than  a-]_.  For  U2^  - Ug  > a-j^,  condition  (20)  can  he  satisfied  only 
with  c'  > 0.  Thus 

0 < c'  < &2}  that  is,  0 < c - Ug  < a2  (22) 

Now,  note  that  condition  ( 20)  is  act-ually  , 

|Ui  - cj  < ai 

and  hence 

. -a^  < UjL  - c < a^  (23) 

Adding  relations  (22)  and  (23),  it  is  found  that 


Ui  - Ug  < + ag 


(21^) 


is  a necessary  condition  for  the  existence  of  suhsonlc  disturhances . 
By  introducing  the  average  velocity  of  sound 


(25) 


and  the  average  Mach  ninnher  of  relative  motion 


M = 


(26) 


defined  with  respect  to  this  average  sound  velocity,  it  may  he  concluded 
that  suhsonic  distxirhances  cannot  exist  (and  the  ng  zone  may  he 
expected  to  he  stable)  if  the  average  Mach  number  M of  the  relative 
motion  exceeds  the  value  2,  as  stated  at  the  heginning  of  this  section. 
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In  terms  of  the  Mach  numbers  of  the  two  streams,  condition  (24-) 
states  that  If 


Ml  > 1 + jr^(i  + Mg) 


(27) 


Stibsonlc  disturhances  cannot  exist,  and  the  motion  Is  stable.  (See 

fig.  3.) 

In  the  case  of  constant  total  enthalpy,  the  temperature  ratio 
T2/T1  caJi  expressed  as 


^ _ 1 + (7  - l)Mi^/2 

l+(r-l)M2^/2 


Then  condition  (27)  becomes 


Ml  >1  + 


1 + (7  - 1)Mi2/2 
1 + (7  - l)J^^/2 


11/2 


(1  + M2) 


The  critical  condition  Is 


(28) 


(29) 


Ml  = 1 + 


1 + (7  - 1)M^2/2 
1 + (7  - 1)M2^/2 


i/2 


(30) 


Removing  the  square  root. 


1 + (7  - 

f2  1 

[Ml  - 1) 

1 + (7  - 1)M22] 

[2  * 

[“2  + 1) 

This  form  suggests  an  obvious  solution,  M2^  = Another  solution 

can  then  be  eeislly  obtained  as 


Ml  = 


(3  - 7)Mg  + k 
(3  - 7)  - 2(7  - 1)M2 


(31) 
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By  subs'tiiTa'bliig  in  equation  (30)^  ^ = -Mg  is  found  to  "be  an  extra- 
neous solution.  The  only  solution  of  equation  (30)  is  then  given  by 
equation  (3I) . 

Thus  the  flow  is  stable  if 


(3  - 7)M2  + h 
M,  > 

(3  - 7)  - 2(7  - 1)>^ 


(32) 


The  curve  for  equation  (3I)  is  shown  in  figure  It,  where  the  region  to 
the  left  of  curve  A is  the  region  of  possible  instability.  Note  the 
symmetry  with  respect  to  the  line  = -Mg,  corresponding  to  a change 

of  the  positive  direction  of  the  x-axis  in  the  physical  problem.  The 
line  = Mg  is  drawn  to  take  care  of  the  condition  U-j_  - Ug  > 0, 

since  U3_  = Ug  when  M^_  = Mg.  The  asymptotes  of  the  curve  for  equa- 
tion (31)  are 


Mg  = (3  - 7)/2(7  - 1) 


and 


Ml  = -(3  •-  7)/2(7  - 1) 


With  7 = l.ll-,  these  values  are  ±2.  The  intercepts  are 


Ml  = V(3  - 7) 


and 


M2  = -V(3  - 7) 


(33) 


{3^) 


(35) 

(36) 


With  7 = 1.4,  these  values  are  +2.5.  (Cf.  reference  9.) 
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GENEE?AL  STUDY  OF  SMALL  DISTUffflANCES  IN  A NEAEU.Y  PAEALLEL 
FLOW  FIELD  IN  A COMPRESSIBLE  FLUID 


The  general  theory  of  a small  dlsturhance  In  a field  of  nearly  par- 
allel flow  of  a gas  has  been  developed  in  reference  1.  The  rigorous 
mathematical  proof  and  improvement  of  the  theory  are  given  in  detail  in 
references  5 and  6.  In  this  report,  merely  the  main  conclusions  and 
their  physical  Interpretations  sire  outlined  without  going  into  the 
details.  Applications  of  the  theory  to  the  specific  case  at  hand  will 
he  discTissed  in  some  detail. 

Consider  a nearly  parallel  stream  with  dimensionless  velocity  and 
ten^ierature  distributions  w(y)  and  T(y) . The  neglect  of  the  depend- 
ence of  these  quantities  on  x and  the  omission  of  the  y-component  of 
the  basic  flow  can  be  justified  by  detailed  investigations.  The  line- 
arized differential  equations  for  small  disturbances  then  possess  solu- 
tions of  the  type 


u' 

v' 

P' 

P' 


Re|f(y)e^(^-'''ti] 

Re 

Re  [^(y)e^(x-ct)^ 

Re  [^(y)e^(^-ct)] 


T' 


Re[0(y)e^(^-'^'tll 


(37) 


where  u',  v',  p’,  p',  and  T'  are  the  perturbations  of  the  two  com- 
ponents of  velocity,  the  density,  the  pressure,  and  the  ten^ierature,  all 
in  a suitably  defined  dimensionless  form.  The  constants  a and  c are, 
respectively,  the  real  wave  number  and  the  complex  wave  speed. 

The  differential  equations  for  the  amplitude  fimctions  f(y),  0(y) ^ 
r(y),  Jr(y)j  aiid  0(y)  are  rather  complicated.  However,  if  the  effects 
of  viscosity  and  heat  conduction  axe  neglected,  they  become  a relatively 
simple  system,  which  can  be  reduced  to  the  following  single  differential 
equation  for  0(y) : 
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(v  - c)0'  - w’0 


d 

<3y|r  - - c)2J 


a^(w  - c)0 


T 


= 0 


(38) 


All  the  other  variahles  can  he  expressed  in  terms  of  0 and  0’  as 
follows : 


f = -i[m^(w  - c)w'0  - T0’  I T - m2(w  - c)^ 
r = ijy0’  + if)  + p'^^(w  - c) 
jt  = l7M^pj^(w  - c)f  + w’^ 


e = T(jt/p  - r/p) 


(39) 


The  boundary  conditions  are  that  the  disturbance  shsuld  be  bounded 
as  y — > ±00, 

The  inviscid  system  would  have  given  a well-defined  characteristic- 
value  problem  if  it  were  not  for  the  fact  that  differential  eq^uation  (38) 
has  a singularity  at  the  point  y = y^,  where  w(y)  = c.  This  singu- 
larity disappears  only  when  d(pw;')/ly  = 0 at  the  same  point.  Other- 
wise, a solution  of  equation  (38)  has  a logarithmic  singularity  at 
y = y^  and  the  characteristic-value  problem  associated  with  this  equa- 
tion becomes  indeterminate  until  the  proper  branch  of  the  solution  is 
determined . 

The  determination  of  the  proper  branch  of  the  solution  and  its 
associated  physical  interpretation  is  one  of  the  most  delicate  points 
in  the  theory  of  hydrodynamic  stability.  The  mathematical  analysis  of 
the  solutions  of  the  complete  viscous  equations  and  their  limiting  solu- 
tions will  be  made  first  before  discussing  their  physical  Interpretation. 

The  complete  system  of  viscous  equations  can  be  shown  to  be  equiva- 
lent to  a system  of  six  linear  equations  of  the  first  order  in  six 
unknowns.  Thiis,  there  are  six  independent  solutions.  These  solutions 
have  been  formally  obtained  as  asyBg)totlc  series  in  reference  1,  and 
their  rigorous  mathematical  investigation  has  been  carried  out  in 
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references  5 and  6.  It  is  foTond  that  two  of  the  six  solutions  can  he 
expressed  in  asynptotic  series  of  the  form 


f = f(®)  + — f(^)  + . . . 


0 = 0(0)  + — 0(1)  + . . , 


(to) 


where  f(0)^  0^^^^  • • • are  the  inviscid  solutions  satisfying  equa- 

tions (38)  and  C39)  and  X^  = oR,  where  R is  the  Reynolds  number  based 
on  the  thickness  of  the  wiixing  region.  Thus,  the  formal  limit  of  equa- 
tions (ho)  does  approach  the  inviscid  solution,  hut  a complete  study  of 
these  equations  also  carries  the  knowledge  of  the  proper  branch  to  he 
■used. 


Four  other  solutions  of  the  complete  system  of  viscous  eq.uations 
are  of  the  form 


where 


f = F exp^XQj^^ 
0=0  exp(XQ^) 


i = 1,  2,  3,  4 ^ 


(to) 


Qi  = -Qg 


(w  - c)  dy 


(to) 


and  F and  0 can  he  expressed  as  power  series  of  X,  involving  only 
a finite  number  of  positive  powers. 
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In  the  establishment  of  these  asymptotic  solutions,  it  is  shown  that 
the  lines  Re^Q^^  = 0,  Re^Q^^  = 0,  and  Re^Q^^  - = 0 are  of  interest. 

The  geometry  of  these  lines  relative  to  the  point  (where  w = c)  and 

the  real  axis  of  the  y-plane  are  shown  in  figure  5.  There  are  asymptotic 
solutions  expressed  in  equations  (l|-0)  and  (^4-1)  which  maintain  the  same 
analytical  expression  on  the  two  sides  of  the  dotted  lines.  However,  in 
crossing  the  solid  lines,  they  generally  change  their  behavior.  Thus, 
the  following  conclusion  may  be  drawn:  The  proper  branch  of  the  multiple- 

valued  asyn5)totic  solutions  is  obtained  by  taking  a path  in  the  con5)lex 
y-plane  below  the  point  y = y^  (in  this  case  w(y)  is  monotonically 

increasing  along  the  real  axis) . 

This  is  the  breinch  taken  in  reference  1.  By  examining  the  behavior 
of  solutions  of  the  type  given  in  equations  (i|-l),  it  can  be  easily  shown 
that  they  diverge  for  either  positive  or  negative  large  values  of  y. 

These  solutions  should,  therefore,  be  rejected  in  the  present  problem. 

The  effect  of  viscosity  is  then  to  be  obtained  through  solutions  of  the 
form  of  eqmtions  (40) . 


GENERAL  STABILITY  CHARACTERISTICS  IN 
THE  INVISCID  CASE 


Much  of  the  discxission  of  the  inviscid  case  in  reference  1 applies 
to  the  present  case.  However,  as  noted  above,  in  the  present  problem, 
the  possibility  can  be  more  readily  realized  that  subsonic  distvurbances 
may  not  exist  at  all,  and  that  the  motion  may  then  be  e3q>ected  to  be 
completely  stable  with  respect  to  ’small  disturbances.  Another  main 
difference  lies  in  the  "steady"  deviation,  that  is,  solutions  of  equa- 
tion (38)  with  a = 0.  In  the  present  discussion,  the  general  line  of 
discussion  in  reference  1 will  be  followed.  However,  the  difference  in 
the  boundary  conditions  often  catises  a difference  in  the  method  of 
analysis.  The  arguments  are,  therefore,  presented  in  some  detail.  Also, 
it  will  be  shown  that  disturbances  having  wave-lengths  slightly  longer 
than  that  of  the  neutral  subsonic  disturbances  of  finite  wave  length 
are  unstable  while  those  with  slightly  shorter  wave  lengths  are  the 
stable  ones.  This  conclusion  also  applies  to  the  case  of  the  boundary 
layer,  but  it  was  not  obtained  in  reference  1.  The  analysis  also  leads 
to  an  approximate  estimation  of  the  dependence  of  eimplification  on  wave 
length.  This  will  be  used  for  the  calcizlation  of  the  an5>lification  of 
the  disturbances  in  the  section  "Self -Excited  Oscillations." 

Some  general  analytical  properties  of  the  solutions  of  equation  (38) 
will  be  first  sximmarlzed,  particularly  for  the  case  of  subsonic  dist-urb- 
ances.  In  this  case. 
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T - - c)^  > 0 (43) 

for  both  free  streams.  For  large  positive  values  of  j,  equation  (38) 
may  he  approximated  hy 


0'  ' - = 0 with  1 - - c)^ 


(44) 


For  large  negative  values  of  y,  it  may  he  approximated  hy 


0"  - p2^0  = 0 


a 


(45) 


Since  hoth  and  ^2^  positive,  the  solution  0(y)  is  e3?ponen- 

tial  in  nature  for  large  values  of  y in  the  case  of  subsonic  disturb- 
ances, with  the  exception  of  the  case  a = 0.  In  that  case,  t^ro  inde- 
pendent solutions  of  equation  (38)  are 


<Ii2_  = w - c 


(46) 


and 


<I) 


2 


(w  - 


(47) 


The  first  solution  is  hounded  while  the  second  varies  linearly  with  y 
for  large  valiies  of  y. 

For  any  value  of  c,  the  hounded  solution  (equation  (46))  corre- 
sponds to  no  disturbance  at  infinity.  In  fact,  for  a = 0,  the  disturb- 
ance v'  is  identically  zero,  hy  equation  (37)-  The  other  components 
of  the  disturbance  are  given  hy  equations  (39)^  and  it  can  he  easily 
verified  that  they  all  vanish  at  large  distances. 

It  can  he  shown  from  the  general  nature  of  the  temperature  and 
velocity  distributions  that  condition  (43)  holds  throughout  the  mixing 
zone  if  it  holds  in  hoth  free  streams.  It  is  then  obvious  that 
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eq-uatlon  (38)  has  a singularity  only  at  w - c = 0.  Actually,  even 
if  T - - c)^  = 0,  it  only  gives  rise  to  an  apparent  singularity. 

Analysis  of  the  solution  in  the  neighbor  hood  of  the  point  y = y^, 
(where  w = c)  gives  the  following  two  solutions: 


01  = (y  - yc)si(y  - yc) 


02  = S2(y  - yc)  + ^01  (y  - yc) 


where 


(48) 


(49) 


fmfi  g^  and  g2  are  power  series  in  (y  - y^,)  with  g]_(0)  = 0 0 

flri(^  63  = ^c/^c'  0 proper  branch  of  tbe  logarithmic  function 

is  to  be  taken  in  accordance  with  the  method  discussed  in  the  preceding 
section. 

For  . real  values  of  c,  it  can  be  shown  that  the  Reynolds  shear 
stress  ' 


T = -pu'v’ 


(50) 


is  a constant  except  for  a possible  jump  at  y = y^-  In  fact,  this 
jump  is 

[T]  (51) 

when  one  passes  from  - 0 to  y^.  + 0.  Thus,  for  a neutral  disturb- 
ance with  a 0,  the  condition 

K = 0 (52) 


must  be  satisfied,  since  t = 0 for  y — >-  ±i»  and,  therefore,  does  not 
have  any  junp. 
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Thus,  for  the  existence  of  a sxibsonic  dlsturhance,  the  quantity. 


(53) 


must  vanish  at  some  point  in  the  field;  furthermore,  the  corresponding 
value  Vg  of  w must  he  subsonic  relative  to  both  streams;  that  is. 


1 - 


M. 


< Wo 


(5k) 


The  above  reasoning  cein  be  applied  to  the  case  of  the  boundary  layer; 
a somewhat  different  argument  wan  used  in  reference  1.  A little"  calcu- 
lation will  show  that,  in  general,  condition  (53)  will  be  satisfied, 
although  condition  (54)  may  not.  If  the  latter  is  also  satisfied,  then 
there  actually  exists  a subsonic  disturbance  with  c = Wg.  To  prove 
this,  eqmtion  (38)  is  rewritten  in  the  form 


where 


(55) 


ii(y)  = [t  - "^  > 0 


(56) 


With  c = Wg,  q.(y)  is  regular  all  along  the  real  axis.  This  is  a 

P 

characteristic-value  problem  for  the  parameter  k = -or  and  is  associ- 
ated with  the  variational  principle 


dy  = 0 


(57) 
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with 


^ dy  = Constant 


(58) 


The  least  value  for  -a^  is  then  given  hy  the  minimum  value  of  the 
ratio 


(59) 


for  1 functions  f (y)  such  that  the  integrals  in  this  ratio  are  con- 
vergent. Now,  it  is  necessary  only  to  show  that  for  certain  fimctions 
the  ratio  of  the  integrals  is  negative. 

In  the  choice  of  such  a function,  it  is  necessary  for  only  the 
numerator  in  equation  (59)  to  he  convergent.  The  convergence  of  the 
integral  in  the  denominator  is  immaterial,  for  one  can  always  modify 
the  function  at  sufficiently  large  values  of  ly]  so  that  the  denomi- 
nator becomes  convergent  without  altering  the  sign  of  the  numerator. 

To  choose  a function  so  that  the  integral 


roo 

(hf '2  + af^)  dy  (60) 


is  negative,  it  is  first  noted  that 


1 = 0 for  f = w - c (61) 

This  follows  from  the  fact  that  0 = w - c is  a solution  of  equation  (55) 
when  a = 0.  It  can  also  he  directly  verified.  Obviously,  the  value  of 
I is  not  changed  if  f is  now  changed  to  Iw  - cl . A further  modifi- 
cation of  the  function  wo\ald  yield  the  desired  resTilt.  In  the  neighbor- 
hood of  yc,  1 w - c I is  small,  but  w'  is  finite.  If  in  this  neighbor- 
hood the  function  f(y)  = [w  - c)  is  replaced  by  a horizontal  straight 
line,  the  integrand  in  equation  (60)  is  certainly  decreased  (since  | >0 
pnd  a change  in  f ' causes  a larger  change  in  the  Integral  than  that 
in  f),  and  the  resultant  integral  is  negative.  This  completes  the  proof 
desired. 
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Thus,  the  sufficiency  of  conditions  (53)  eind  (54)  for  the  existence 
of  subsonic  disturbances  is  also  established.  Let  the  corresponding 
value  of  a be  denoted  by  Og. 

Rext,  it  will  be  shown  that  for  a slightly  less  than  Og  the 
disturbances  are  unstable.  To  prove  this  the  quantity  (dc/dk)g  will 
be  calculated,  a-nfl  its  imaginary  part  will  be  shown  to  be  positive. 

Consider  a characteristic  function  0(y;k,c)  of  equation  (55) 
corresponding  to  a given  value  of  k-  As  c changes,  k also  changes. 
If  eq.uation  (55)  is  differentiated  with  respect  to  k,  the  following 
eq.uation  is  obtained 


3-y 


^ + ^ = 
dk  T 


0 


(62) 


where 

+ (63) 

8k  8c  dk 


Row,  multiply  equation  (62)  by  0,  eq,mtion  (55)  by  sxibtract  the 
results,  and  then  integrate  with  respect  to  y along  a path  in  the  com- 
plex plane  which  leads  from  y = -00  on  the  real  axis  to  y = «>  on  the 
real  axis  but  passes  below  the  point  y^..  The  following  equation  is  then 
obtained  for  dc/ dk: 


^ r ^ 02 

dk  J_„  ^c\dy/  8c  _ 

So  far  equation  (64)  is  general.  Row,  specialize  to  the  case  of 
the  neutral  distiu’bance  in  question.  Most  of  the  integration  can  then 
be  carried  out  along  the  real  axis,  with  real  res-ultant  values.  However, 
the  integral 


dy 


=/: 


K 

T 


(64) 
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must  "be  evaluated  along  an  indented  path  vith  a small  detour  "below 

since  8q/8c  has  a pole  at  that  point.  Calculations  show  that  the 
imaginary  part  of  J is 


Im(j) 


(66) 


if  this  quantity  does  not  vanish.  Indeed^  this  quantity  can  he  seen 
to  he  negative.  Thiis^  equation  (6l|-)  yields  a relation  of  the  form 


— (P  + IQ)  = E > 0 

dk 


where  Q = Im(j) . Thus^ 


dc 

dk 


E 

p2  + q2 


(P  - iQ) 


has  a positive  imaginary  part,  as  required. 

In  the  incompressible  case,  equation  (64)  reduces  to 


(67) 


dc 

dh 


(w  - c)2 


02  dy 


and  the  imaginary  part  of  dk/dc  is 


(68) 


Im( dk/dc)  = Jt^0c|’^c')^  "^c"'  f 


(69) 
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WEUTRAL  OSCniATIONS 


From  the  discussions  of  the  last  section,  another  limitation  is 
ingjosed  on  the  occurrence  of  neutral  subsonic  disturbances.  Since  the 
wave  speed  of  the  neutral  disturbance  must  be  equal  to  the  flow  speed 
Wg  at  the  point  where 


d / dw\ 
dy  \ dy/ 


= 0 


(70) 


subsonic  disturbeinces  cannot  occur  if  the  speed  Wg  is  supersonic  rela- 
tive to  either  stream.  Calculations  of  Wg  are  made  for  the  case  of 
uniform  total  enthalpy.  Instead  of  using  equation  (70)  directly,  it 
is  fovmd  convenient  to  transform  it  into  the  form 


-1 

= 0 with  Tjr  = f(^)  and  w = f'(^)  (71) 


where  -ijf  is  the  stream  function,  since  the  relation  w(ijf)  is  identical 
with  that  in  the  incompressible  case.  The  results  of  these  calculations 
aore  shown  in  table  m.  In  con^iaring  these  values  against  condition  (54), 
it  is  found  that  the  condition 


1 + 


(7  - l)Mi" 


- w2 


1 

/ 


M-, 


is  never  violated,  but  the  values  of  Wg  below  the  solid  horizontal 
lines  in  the  table  are  too  high  to  satisfy  the  condition 


Wc 


Thus,  for  a given  Ifeich  mmiber  in  the  stream  Mg,  the  Mach  number 
in  the  other  stream  can  become  so  high  that  the  speed  of  flow  corre- 
sponding to  condition  (70)  becomes  supersonic  relative  to  the  slower 
free  stream.  Subsonic  dlstijrbances  then  cannot  exist. 
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The  critical  case  of  sonic  disturhance  is  reached  when 


Calculations  for  this  case  are  tabulated  in  table  IV.  The  corresponding 
values  of  Mp  and  M2  are  plotted  as  curve  B in  figttre  k to  mark  the 
limit  of  stability.  Only  the  region  to  the  left  of  this  curve  can  have 
subsonic  disturbances.  It  is  seen  that  the  condition  is  more  restrictive 
than  that  obtained  from  general  considerations  alone.  In  particular,  in 
the  case  where  one  stream  is  at  rest,  the  cxirve  B shows  that  the  flow 
becomes  con^jletely  stable  at  a free-stream  Mach  number  of  l.Tj  ia  con- 
trast to  the  value  2.5  given  by  curve  A.  Since  they  are,  however,  not 
very  much  different,  it  may  be  simnised  that  the  exact  distributions' of 
ten^jeratiire  and  velocity  in  the  mixing  region  may  not  be  too  in^jortant 
in  determining  the  over-all  stability  characteristics. 

It  may  be  recalled  that  the  restriction  to  the  case  of  constant 
total  enthalpy  can  be  immediately  removed  by  considering  a moving  obser- 
ver. Thus,  for  all  cases  with  Prandtl  number  equal  to  unity,  it  is  nec- 
essary only  to  convert  the  values  of  the  speeds  Up,  U2,  and  Wg 

Involved.  Condition  (70)  is  not  modified  by  the  reference  to  a moving 
observer. 

With  the  neutral  wave  speed  thus  determined,  eofuation  (55)  can  be 
integrated  to  give  the  amplitude  of  the  oscillations.  For  this  purpose, 
it  is  necessary  to  find  the  proper  value  for  a.  This  can  be  done  by 
several  trials,  with  the  first  approximation  given  by  the  ratio  of  the 
integrals  in  eq-uation  (59) . The  ratio  will  yield  the  characteristic 
value  only  when  the  function  f (y)  is  the  characteristic  function,  but 
it  is  known  that  any  reasonable  approximation  to  it  will  give  a very 
close  approximation  to  the  characteristic  valtie. 

Calculations  of  these  neutral  oscl2JLations  are  carried  out  for 
several  Mach  numbers  in  one  stream  with  the  other  at  rest.  These  are 
given  in  table  Y and  figure  6. 


SELF-EXCITED  OSCILLATIONS 


The  formulas  given  in  the  section  "General  Stability  Characteristics 
in  the  Invlscld  Case"  have  been  used  to  calculate  the  characteristics  of 
self-excited  oscillations.  It  is  found  that 


— = 0.177  - 0.2091 
da 


2k 
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for  the  case  M]_  = 1 eind 

— = 0.093  - 0.2871 

do, 

for  the  incompressihle  case.  The  results  are  shown  dieigraiiimatically  in 
figvire  7.  It  is  seen  that  the  extent  of  anipllfication  Is  fairly  large. 
There  is  also  indicated  a decrease  of  an5>lification  with  increasing  Mach 
number. 

It  would  he  easy  to  calculate  the  amplitude  function  of  these  self- 
excited  oscillations  hy  using  the  characteristic  values  obtained  above. 
This  was  not  carried  out  becaoise  of  limitations  in  time. 


NATURE  OF  OSCILLATIONS  IN  LIMIT  OF 
INFINITE  REYNOLDS  NUMBER 


Calculations  in  the  section  "Neutral  Oscillations"  Indicate  that 
the  neutral  oscll lation  has  quite  a singjle  amplitude  distribution.  In 
fact,  it  does  not  show  any  node.  The  amplified  distmrbances  are  expected 
to  show  similar  characteristics.  The  mathematical  theory  (references  5 
and  6) , however,  indicates  that  the  danped  oscillations  behave  in  a much 
more  conplicated  manner.  It  is  concluded  that  there  is  always  a finite 
viscous  region  in  the  interior  of  the  fluid,  no  matter  how  large  the 
Reynolds  number  may  be.  In  fact,  a minimum  width  of  this  region  is 
determined.  This  has  to  do  with  the  crossing  of  the  solid  lines  by  the 
real  axis  of  y in  figure  5.  The  solution  in  the  part  of  the  real  eixis 
between  the  solid  lines  shows  exponential  behavior  - and  therefore  vIbcoub 
nature  - if  the  solution  in  the  outside  parts  shows  the  invLscid  behavior. 
This  type  of  conclusion  has  been  reached  in  reference  1.  However,  it 
was  possible  only  to  s-uggest  that  such  viscous  behavior  would  occur  at 
the  solid  lines.  The  inproved  theory  shows  that  it  must  occur  through- 
out the  region  in  between.  It  is  Important  to  note  tlat  the  complex  con- 
jugate of  an  amplified  solution  does  not  represent  a damped  oscillation 
and  vice  versa,  althotigh  this  conclusion  can  be  easily  reached  by  a 
cursory  examination  of  the  invlscid  equations.  The  damped  oscillations 
do  not  satisfy  the  invlscid  equation  a1 1 along  the  real  axis;  otherwise, 
they  could  not  tahe  on  the  proper  branch  of  the  logarithm  as  specified 
in  the  section  "General  Study  of  Small  Dlstinrbances  in  a Nearly  Peirallel 
Flow  Field  in  a Conpresslble  Fluid."  They  exhibit  a behavior  very  much 
like  that  of  the  vorticity  distribution  in  fully  developed  turbulence 
flow.  For  large  Reynolds  numbers,  there  is  one  part  of  the  space  where 
the  vorticity  is  highly  concentrated;  in  another  part,  there  is  very 
little  vorticity.  This  1 1.1  \istrates  the  two  kinds  of  limiting  behavior 
of  a viscous  fluid  in  the  limit  of  infinite  Reynolds  mmiber:  In  one 

paxt  of  the  field,  the  invlscid  behavior  is  approached;  in  another  part, 
it  becomes  highly  oscillatoiy  spatially. 
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CONCLUSIONS 


From  a study  of  the  stahility  of  the  mlxiriig  of  two  parallel  streams 
In,  a gas^  the  following  conclvisions  may  he  drawn: 

1.  If  the  relative  speed  of  the  two  parallel  streams  exceeds  the  sum 
of  their  velocities  of  sound,  subsonic  oscillations  cannot  occxor,  and 
the  mixing  region  may  be  expected  to  be  stable  with  respect  to  small 
disturbances . 

2.  A further  necessary  condition  for  the  possible  occurrence  of 
small  subsonic  dlstirrbeinces  is  that  somewhere  in  the  field 


where  y is  a posit ionad.  coordinate  across  the  stream,  p is  the  den- 
sity of  the  gas,  and  w is  a dimensionless  velocity  distribution.  This 
condition  is  usually  satisfied  for  the  present  class  of  problems. 

3.  If  the  speed  of  the  flow  at  such  a point  is  denoted  by  Wg, 
then  the  field  of  flow  can  execute  a neutral  wavy  oscillation  having 
a finite  wave  length  and  propagating  with  the  speed  c = Wg  if  and 
only  if  Wg  is  subsonic  relative  to  both  streams.  There  is  no  other 

possible  neutral  oscillation.  This  leads  to  a more  strict  condition 
of  stability  than  that  given  by  conclusion  1. 

4.  Under  the  above  conditions,  the  field  of  flow  can  execute 
amplified  wavy  oscillations  having  wave  lengths  longer  than  that  of 
the  neutral  oscillation.  Oscillations  having  shorter  wave  lengths  axe 
damped.  (This  specific  form  of  the  conclusion  was  not  obtained  for 
the  boundary  layer  at  a solid  surface  in  NACA  TN  1 1 15,  but  its  validity 
can  be  shown  by  the  present  method.)  The  extent  of  arapl±ficatlon  in 
such  cases  is  fairly  large. 

5.  At  large  Reynolds  numbers,  the  amplified  disturbances  are 
essentially  free  from  the  effect  of  viscosity.  On  the  other  hand, 
disturbances  with  finite  damping  are  expected  to  exhibit  a highly 
oscillatory  behavior  over  a finite  region  in  the  field  of  flow.  This 
is  similar  to  the  structure  of  the  vortlcity  field  in  fully  developed 
turbulence . 

6.  For  the  case  of  constant  enthalpy  with  one  stream  at  rest,  the 
wave  length  of  the  neutral  disturbances  increases  with  increasing  Mach 
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number  in  ■fclie  other  stream.  The  flow  becomes  completely  stable  at  a 
free-stream  Mach  nxmiber  of  1.7.  This  is  more  restrictive  than  the 
value  2.5  obtained  by  applying  conclusion  1 to  the  present  case. 


Massachusetts  Institute  of  Technology 
Cambridge,  Mass.,  July  30,  1952 
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APPENDIX 

SYMBOLS 


The  quantities  'bearing  subscript  1 in  the  last  column  are  the 
dimensional  q.uantitles  in  the  first  free  stream.  Corresponding  quan- 
tities in  the  second  stream  hear  a subscript  2.  The  quantities  with- 
out a prime  satisfy  the  equations  of  steady  motion;  those  with  primes 
satisfy  the  disturbance  eqmtlons. 

Dimensional  Dimensionless  Reference 

q.uantlties  quantities  quantities 

Positional  cooi-dinates ; 

X*  X I = 2\2ViX*/Ui 


y* 

y = rj2.^ 

1 

Time; 

t* 

t 

Z/Ui 

Velocity  components 

in  directions  of  x-  and  y-axes. 

respectively: 

u*  + u-^' 

w(y)+f(y)e^(^-‘''t) 

Ul 

Y*  ^ Y*  * 

( y ) ® ^ 

Ui 

Density  of  gas; 

p*  + p*' 

p(y)+r(y)e^(=‘-'=*) 

Pi 

Pressinre  of  gas; 

p*  + p*’ 

p(y)+rt(y)ei^(x-'='‘^) 

Pi 

Tenqperattire  of  gas; 

T*  ^ Qi*t 

T(y)+0(y)ei^(^-''*) 

Tl 

Coefficient  of  viscosity  of  gas; 
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DlmenslonEil  Dimfinslonless 

quemtities  quantities 

Thermal  conductivity: 

k*  p/Pr 


Wave  number  of  dlsturhance: 

a*  = 2ir/X*  a = 2jt/X 

Phase  velocity  of  dlsturhance: 

c*  c 


Specific  heat  at  constant  volume: 


^eclflc  heat  at  constant  pressure:- 

cp  r 


Gas  consteint  per  gram: 

R*  7-1 

Reynolds  number 

R = PiUil/wi 

Mach  number 

Ml  = Ui|j7H*Ti 


Reference 

quantities 


1-1 


Ul 


Cv 


Cv 


Prandtl  number 


Pr  = Cp(i*/k*^ 
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TABLE  I 

VELOCITy  DISTRIBUTIONS  IN  INCOMPRESSIBLE  CASE  - Continued 


(d)  For  X = 0.8  and  6 = 0.860 


V 

n 

w 

n 

w 

-9-0 

0.111 

-2.6 

0.209 

1.4 

0.813 

-8.8 

.111 

-2.4 

.224 

1.6 

.842 

-8.4 

.111 

-2.2 

.240 

1.8 

.869 

-8.0 

.111 

-2.0 

.258 

2.0 

.893 

-7.6 

.112 

-1.8 

.279 

2.2 

.9lli- 

-7.2 

.112 

-1.6 

.301 

2.4 

.932 

-6.8 

.113 

-1.4 

.326 

2.6 

.947 

-6.4 

.115 

-1.2 

.352 

2.8 

.960 

-6.0 

.117 

-1.0 

.382 

3.0 

.970 

-5.6 

.120 

-.8 

.413 

3.2 

.977 

-5.2 

.124 

-.6 

.446 

3A 

.984 

-4.8 

.129 

-.4 

.481 

3.6 

.988 

-4.4 

.137 

-.2 

.518 

3.8 

.992 

-4.0 

.146 

0 

.556 

4.0 

.994 

-3.8 

.152 

.2 

.594 

4.2 

.996 

-3.6 

.159 

.4 

•633 

4.4 

.998 

-3.^ 

.167 

.6 

.672 

4.6 

.999 

-3.2 

.175 

.8 

.709 

4.8 

.999 

-3.0 

.185 

1.0 

.746 

5.0 

1.000 

-2.8 

.196 

1.2 

.780 
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TABLE  II 

VELOCITY  AND  TEMPERATURE  DISTRIBUTIONS 
IN  COMPRESSIBLE  CASE 


[toe  stream  at  rest;  X = 1;  0 


(a)  = 0.5 


■H 

w 

T 

n 

w 

T 

-12.6 

-12.2 

-11.8 

0 

1.050 

-1.4 

-1.2 

-1.0 

0.261 

.289 

.320 

!-■  H H 

bob 

-11.4 

-n.o 

-10.6 

-.8 

-.6 

-.4 

.352 

.386 

.423 

1.044 

1.043 

1.041 

-10.2 

-.2 

.461 

1.039 

-9.8 

0 

.500 

1.038 

-9A 

0 

1.050 

.2 

.540 

1.035 

-9.0 

.001 

.4 

.581 

1.033 

-8.6 

-8  2 

.001 

.002 

.003 

.6 

.8 

1.0 

.622 

.662 

.702 

1.031 

1.028 

1.025 

-7.8 

-7.1^ 

.004 

1.2 

.740 

1.023 

-7.0 

.006 

1.4 

.776 

1.020 

-6.6 

.009 

1.6 

.809 

1.017 

-6.2 

.013 

1.8 

.840 

1.015 

-5.8 

.019 

2.0 

.868 

1.012 

-5.k 

.025 

2.2 

.892 

1.010 

-5.0 

.033 

1.049 

2.4 

.914 

1.008 

-4.6 

.043 

2.6 

.932 

1.007 

-4.2 

.054 

2.8 

.947 

1.005 

-3.8 

.069 

3.0 

.960 

1.004 

-3.4 

.086 

3.i^ 

.978 

1.002 

-3.0 

.108 

3.8 

.989 

1.001 

-2.8 

.121 

4.2 

.995 

1.001 

-2.6 

.136 

4.6 

.998 

1.000 

-2.4 

.152 

5.0 

.999 

1.000 

-2.2 

.170 

5.4 

1.000 

1.000 

-2.0 

.190 

1.048 

5.6 

1.000 

1.000 

■a 

.211 

.235 

1.048 

1.047 

6.0 

1.000 

1.000 

w 


T 


V 


T 


w 

T 

w 

0 

1.200 

^9 

0.282 

.309 

.336 

.366 

.397 

.430 

0 

1.200 

-.2 

.465 

.001 

0 

.500 

.002 

.2 

.537 

.002 

.4 

.574 

.004 

.6 

.611 

.005 

.8 

.648 

.007 

1.0 

.685 

.010 

1.2 

.721 

.014 

1.4 

.756 

.019 

1.6 

.789 

.024 

1.8 

.820 

.031 

2.0 

.848 

.039 

2.2 

CO 

• 

.048 

2.4 

.897 

.059 

2.6 

.918 

.072 

1.199 

2.8 

.935 

.088 

1.198 

3.0 

.950 

.108 

1.197 

■ 3-^ 

.971 

.132 

1.197 

3.8 

.985 

.145 

1.196 

4.2 

.993 

.160 

1.195 

4.6 

.997 

.177 

1.194 

.999 

.195 

1.192 

1.000 

.214 

1.191 

.235 

.258 

1.189 

1.187 

6.0 

l.lSlt 

1.181 

1.177 

1.173 

1.168  ' 

1.163 

1.157 

1.150 
1.142  ■ 

1.134 

1.125 

1.116 

1.106 

1.096 

1.086 

1.076 

1.066 

1.056 

1.047 

1.039 

1.032 

1.025 

1.020 

1.011 

1.006 

1.003 

1.001 

1.001 

1.000 
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TABLE  II 

VELOCITY  AMD  TEMPERATURE  DISTRIBUTIOMS 


IN  COMPRESSIBLE  CASE  - Concluded 


(c)  Mjl  = 1.5 


BH 

V 

T 

4 

w 

T 

-12.6 

0 

l.ltBO 

-1.4 

0.310 

1.407 

-12.2 

.001 

-1.2 

• 334 

1.400 

-11.8 

.001 

-1.0 

.359 

1.392 

-u.lt- 

.002 

-.8 

.385 

1.384. 

-11.0 

.002 

-.6 

.412 

1.374 

-10.6 

.003 

-.4 

.440 

1.363 

-10.2 

.OOlt- 

-.2 

.470 

1.351 

-9.8 

.006 

0 

.500 

1.338. 

-9.h 

.007 

.2 

.531 

1.323 

-9.0 

.009 

.4 

.564 

1.307 

-8.6 

.012 

.6 

.596 

1.290 

-8.2 

.016 

— 

.8 

.629 

1.272 

-7.8 

.020 

1.0 

.662 

1.253 

-7.1^ 

.025  • 

1.2 

.694 

1.233 

-7.0 

.030 

1.4 

.727 

1.213 

-6.6 

.037 

l.¥^9 

1.6 

.758 

1.192 

-6.2 

.0lt4 

l.kh9 

- 1.8 

.788 

1.171 

-5.8 

.052 

l.4it-9  ■ 

2.0 

00 

• 

1.150 

-5.k 

.062 

l.W 

2.2 

.844 

1.130 

-5.0 

• 073 

I.IA8 

2.4 

.869 

1.111 

-k.6 

.086 

l.itli-7 

2.6 

.891 

1.093 

-It-. 2 

.102 

l,kh5 

2.8 

.911 

1.076 

-3.8 

.120 

l.Wt- 

3.0 

.929 

1.062 

-3-h 

.llt-2 

l.i^H 

3.4 

.957 

1.038 

-3.0 

.167 

l.lt38 

3.8 

.976 

1.021 

-2.8 

.180 

1.435 

4.2 

.988 

1.011 

-2.6 

.196. 

1.433 

4.6 

.994 

1.005 

-2.lt- 

.212 

1.430 

5.0 

.998 

1.002 

-2.2 

.229 

1.426 

5.4 

.999 

1.001 

-2.0 

.2hQ 

1.422 

5.6 

1.000 

1.000 

-1.8 

.267 

1.418 

6.0 

1.000 

1.000 

-1.6 

.288 

1.413 
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TABLE  HI 

WAVE  SPEEDS  OF  POSSIBLE  HEUTRAL  SUBSONIC  DISTURBANCES 


u*/U3^  defined  By  equation  (70^ 


X 

0.2 

0.4 

0.6 

0.8 

1.0 

0 

0.83l». 

0.732 

0.657 

0.613 

0.576 

.5 

.81^5 

.743 

.678 

.632 

.594 

1 

.851 

.762 

.704 

.664 

.633 

2 

.878 

.823 

.796 

.772 

.762 

5 

.9h9 

.940 

.941 

.938 

.938 

10 

.984 

.981 

.981 

.982 

.982 

TABLE  IV 

CONDITIONS  ASSOCIATED  WITH  NEUTRAL  SONIC  DISTURBANCE 


Conditions 

X 

0.2 

0.4 

0.6 

0.8 

1.0 

Ml 

>10 

4.0 

2.7 

2.0 

1.7 

Mg 

>1.92 

.90 

.44 

.167 

0 

c 

sjl 

.91 

.82 

.78 

.73 

NACA  TN  2887 


39 


TABLE  V 

WAVE  NUMBER  AND  AMPLITUDE  OF  NEUTRAL  OSCILLATIONS 

Qr  = e = 1.160J 

(a)  Wave  number  of  neutral  oscillations 
for  various  Mach  numhers 


% 

Wave  number,  o0 

0 

0.i^59 

.5 

1.0 

.371^ 

1.5 

.324 
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TABLE  V 

WAVE  MJMBER  AM)  AMPLITUDE  OF  NEUTRAL  OSCILLATIONS  - Continued 


(L)  Aniplitude  of  neutral  oscillations 
for  = 0 


JK 
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TABLE  V 

WAVE  NUMBER  AND  AMPLITUDE  OF  NEUTRAL  OSCILLATIONS  - Concliided 

(c)  Amplitude  of  neutral  oscillations 
for  = 1 


y 

0 

d0/dy 

y 

0 

d0/dy 

y 

0 

d0/dy 

-l.k 

0.003 

0.003 

-0.3 

0.551 

0.284 

2.9 

0.122 

-0.106 

-7.0 

.005 

.003 

-.2 

.579 

.273 

3.0 

.112 

-.098 

-6.6 

.006 

.005 

-.1 

.605 

.254 

3-1 

.103 

-.090 

-6.2 

.008 

.006 

0 

.629 

.227 

3.2 

.094 

-.082 

-5.8 

.011 

.008 

.1 

.650 

.190 

3.3 

.087 

-.076 

-5.4 

.015 

.011 

.2 

.667 

.142 

3.4 

.079 

. -.070 

-5.0 

.020 

.015 

.3 

.678 

.085 

3.6 

.067 

-.059 

-4.6 

.026 

.020 

.4 

.684 

.020 

3.8 

.056 

-.049 

-4.2 

.035 

.026 

.5 

.682 

-.050 

4.0 

.048 

-.042 

-3.8 

.048 

.035 

.6 

.674 

-.121 

4.2 

.040 

. -.035 

.064 

.047 

.7 

.658 

-.189 

4.4 

•.034 

-.030 

-3.0 

.086 

.063 

.8 

.636 

-.249 

4.6 

.029 

-.025 

-2.6 

.115 

.084 

.9 

.609 

-.297 

4.8 

.024 

-.021 

-2.2 

.153 

.112 

1.0 

.577 

-.331 

5.0 

.020 

-.018 

-2.0 

.177- 

.129 

1.1 

.543 

-.351 

5.2 

.017 

-.015 

-1.9 

.191 

.138 

1.2 

.507 

-.359 

5.4 

.014 

-.013 

-1.8 

.205 

.148 

1.3 

.471 

-.356 

5.6 

.012 

-.011 

-l.f 

.220 

.158 

1.4 

.436 

5.8 

.010 

-.009 

-1.6 

.237 

.169 

1.5 

.403 

-.328 

6.0 

.009 

-.008 

-1.5 

.254 

.180 

1.6 

.371 

-.308 

6.2 

.007 

-.006 

-1.4 

.273 

.192 

1.7 

.341 

-.288 

-1.3 

.292 

.204 

1.8 

•3-3 

-.267 

-1.2 

.313 

.217 

1.9 

.288 

-.246 

-1.1 

.336 

.229 

2.0 

.264 

-.227 

-1.0 

.359 

.242 

2.1 

.242 

-.209 

-.9 

.384 

.254 

2.2 

.222 

-.192 

-.8 

.410 

.265 

2.3 

.204 

-.177 

-.7 

.437 

.275 

2.4 

.187 

-.163 

-.6 

.465 

.282 

2.5 

.171 

-.150 

-.5 

.493 

.287 

2.6 

.157 

-.137 

-.4 

.522 

.288 

2.7 

.144 

-.126 

2.8 

.132 

-.116 

(a)  With  both  Btreams  moving.  (b)  WJ,th  one  etream  at  rest. 

Rnri  M2,  Mach  numbers  Ml,  Mach  number  of  moving 

of  tvro  streams.  stream;  M2  = 0;  7 = 1.^. 


Figure  3.-  Stable  zones  for  various  teraperatures  and  Mach  numbers. 
T2/T1,  ratio  of  absolute  temperatures  of  two  streams . 
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(c)  For  Mx  = 1.0 


SJid  aB  = 0.3714-. 


(d)  For  Mx  =■  1.5  and  cg0  = O.32I4-. 
Figure  6.-  Concluded. 
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(■b)  For  Ml  = 1. 

Figure  7‘“  Extenl  of  amplifica'tloii  at  infini'te  Reynolds  number, 

a,  wave  number;  — , wave  length  in  y-unlts;  y = -i—  Tj  = — ; 

“■  2i/2  2\/2  /viX*/Ui 

Rj  Reynolds  number;  c,  dimensionless  con5>lex  wave  speed  + ic^^, 


NACA  - Laijgley  Field,  Vi*. 


